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ABSTRACT 

The gamma-ray energy spectra of bright blazars of the LAT Bright AGN Sam- 
ple (LBAS) are investigated using Fermi-LAT data. Spectral properties (hardness, 
curvature and variability) established using a data set accumulated over 6 months of 
operation are presented and discussed for different blazar classes and subclasses: Flat 
Spectrum Radio Quasars (FSRQs), Low-synchrotron peaked BLLacs (LSP-BLLacs), 
Intermediate-synchrotron peaked BLLacs (ISP-BLLacs) and High-synchrotron peaked 
BLLacs (HSP-BLLacs). The distribution of photon index (F, obtained from a power-law 
fit above 100 MeV) is found to correlate strongly with blazar subclass. The change in 
spectral index from that averaged over the six month observing period is < 0.2-0.3 when 
the flux varies by about an order of magnitude, with a tendency toward harder spectra 
when the flux is brighter for FSRQs and LSP-BLLacs. A strong departure from a single 
power-law spectrum appears to be a common feature for FSRQs. This feature is also 
present for some high-luminosity LSP-BLLacs, and a small number of ISP-BLLacs. It 
is absent in all LBAS HSP-BLLacs. For 3C 454.3 and AO 0235-^164, the two brightest 
FSRQ source and LSP-BLLac source respectively, a broken power law gives the most 
acceptable of power law, broken power law, and curved forms. The consequences of 
these findings are discussed. 



34 
35 
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Introduction 



37 Launched into a low-Earth orbit on June 11, 2008, the Fermi Gamma-Ray Space Telescope 

38 continues providing excellent gamma-ray data for celestial sources. With significant improvement of 

39 sensitivity and bandpass over its predecessors (jAtwood et al.ll2009l ). the main instrument on Fermi 

40 - the Large Area Telescope, or LAT - enables detailed studies of time-resolved broad-band gamma- 

41 ray spectra of a broad ran ge of sources, including active galaxies. As discovered by EGRET on 

42 the Compton Observatory (iHartman et al.lll992l : iFichtel et al.lll994l ). active galactic nuclei (AGN) 

43 showing strong gamma-ray emission are associated with relativistic jets, whose presence was in- 

44 dependently inferred from morphological and variability studies in other bands. The spectra of 

45 such objects in all observable bands are well-d escribed by broa d power-law or curved distributions, 

46 indicating non-thermal emission mechanisms (jBottcherl 120071 ) . Very generally, the overall broad- 

47 band spectral distributions of such jet-dominated AGN, often called blazars have a two-humped 

48 shape, with the low energy (IR-UV) hump attributed to synchrotron emission of energetic electrons 

49 radiating in ma gnetic field and the high energy hump due t o inverse Cqmptq n scattering by the 

50 same electrons (|Ghisellinilll989l : iDermer &: Schlickeiserlll993l : ISikora et al.lll994l ). 
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51 The first list o f such AGN detected by the Fermz-LAT, the LAT Bright AGN Sample (LBAS) 

52 (|Abdo et al.ll2009el ) includes bright, high-galactic latitude > 10°) AGNs detected by the Fermi- 

53 LAT with high significance (Test Statistic TS> 100) during the first three months of scientific 

54 operation. This sample comprises 58 Flat Spectrum Radio Quasars (FSRQs), 42 BLLac-type ob- 

55 jects (BLLacs), two radio galaxies and four quasars of unknown type. This somewhat conventional 

56 classificati on was based on the observed optical emission line equivalent widths and the Ca II break 

57 ratio (e.g., iMarcha et al.lll996l ). Following the mod els used to describe the gamma-ray spectra ob- 

58 tai ned with previous g amma-ray observatories (e.g.. lMattox et al.lll996l ). the early analysis reported 

59 in (|Abdo et al.ll2009d ) was carried out by fitting the gamma-ray spectra at energies above 200 MeV 

60 using a simple power law (PL) model. This analysis revealed a fairly distinct spectral separation 

61 between FSRQs and BLLacs, with FSRQs having significantly softer spectra. The boundary pho- 

62 ton index between the two classes was found to be F ~2.2. It has been suggested (jGhisellini et al 



20091 ) that this separation results from different radiation cooling suffered by the electrons due to 



64 distinct accretion regimes in the two blazar classes. 

65 While adopting such a simple spectral model was sufficient to investigate the source spectral 

66 hardness distribution, a PL model was clearly not the most appropriate choice for some bright 

67 sources which exhibited evident breaks or curvatures in their spectra. The departure of the func 



68 tional form from a PL was investigated in some detail for the bright quasar 3C 454.3 (jAbdo et al 



2009dl ) . which underwent strong activity in the summer of 2008. The change of photon index AF 



70 was observed to be 1.2ib0.2, i.e. greater than the value of 0.5 expected from incomplete cooling of 

71 the emitting electrons. The observed break around 2.2 GeV was ascribed to mirroring a similar 

72 feature in the underlying emitting electron energy distribution; the Klein-Nishina effect was not 

73 ruled out, though the importance of photon-photon pair production requires the gamma-ray emis- 

74 sion region to be close to the supermassive black hole. Clearly, understanding the details of the 

75 spectral break is important for understanding the structure and location of the dissipation region 

76 of jets in active galaxies. 

77 The data first obtained with the EGRET instrument, now refined with Ferr ni, imply that the 



78 high Galactic latitude sky emits quasi-diffus e, uniform gamma-ray background (jSreekumar et al 



19981 : IStrong et al.ll2004al : lAbdo et al.l l2009bl ) . Its isotropy points to its extragalactic nature after 



80 subtraction of a quasi-isotropic gamma-ray emission component from cosmic ray electrons in an 

81 extended galactic halo. Most models account for at least a part of this background as originat- 

82 ing from a large number of unresolved point sources, presumably jet-dominated AGN. Here, the 

83 comparison of spectral properties of various classes of AGN, integrated over their space density 

84 and luminosity, against the integral measurement of the unresolved component should provide ad- 

85 ditional clues regarding their contribution to the extragalactic diffuse background: ca n they make 

86 up th e entire background, is another class of sources, e. g. star-forming galaxies (jFields et al 



20081 ) or is an additional truly diffuse component required? The issue is further complicated by 



88 the apparently different spectral forms of the lumino us AGN associated with quasars as compared 

89 to the lineless BL Lac objects, as already hinted in lAbdo et al.l (2009e). It is thus important to 
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90 determine whether the spectral feature seen in 3C 454.3 is common in all blazars and also whether 

91 it is connected to other blazar properties. 

92 Here we report on the detailed spectral analysis of bright LBAS sources using data accumulated 

93 over the first 6 months of the Fermi-hAT all-sky survey. In Section 2, we present the observations 

94 with the Fermi-LAT; Section 3 briefly discusses the classification scheme used in this paper. Section 

95 4 contains the results regarding the photon index and observed deviations from a pure PL. The 

96 results and their consequences are discussed in Section 5. 



97 2. Observations with the Large Area Telescope 

98 The Fermi-LAT is a pair-conversion gamma-ray telescope sensitive to photon energies greater 

99 than 20 MeV. It is made of a tracker (composed of two sections, front and back, with different 
localization capabilities), a calorimeter, and an anticoincidence system to reject the charged-particle 

01 background. The LAT has a large peak effective area (~ 8000 cm^ for 1 GeV on-axis photons in the 

02 event class "diffuse" considered here), viewing ~ 2.4 sr of the full sky with an angular resolution 
i containment angle) better than ^ 1° at E = 1 GeV (jAtwood et al.ll2009l ). 



The data were collected from 4 Aug. 2008 to 1 Feb. 2009 in survey mode. To minimize 

05 systematics, only photons with energies greater than 100 MeV were considered in this analysis. In 

06 order to avoid contamination from Earth limb gamma-rays, a selection on the zenith angle, < 105°, 

07 was applied. The exposure was constant within 20% for all sources and amounted to about 1.5x10^ 

08 m^s at 1 GeV. 

09 This analysis was performed with the standard analysis tool gtlike, part of the Fermi-LAT 

10 ScienceTools software package (version v9rl2). The first set of instrument response functions (IRFs) 

11 tuned with the flight data, P6_V3_DIFFUSE, was used. In contrast to the preflight IRFs, these 

12 IRFs take into account corrections for pile-up effects. This correction being higher for lower energy 

13 photons, the measured photon index of a given source is about 0.1 higher (i.e. t he spectrum i s softer ) 



14 with this IRF set as compared to the P6_V1_DIFFUSE one used previously in lAbdo et al.l ()2009d ). 

15 Photons were selected in circular regions of interest (ROI), 7° in radius, centered at the positions 

16 of the sources of interest. The isotropic background (the sum of residual instrumental background 

17 and extragalactic diffuse gamrna-ray background) was modeled with a simple power-law. The 

18 GALPROP model ( Strong et al. 2004al fbl). version "glLiem_v01.flt", was used for the galactic diffuse 



19 emission, with both flux and spectral photon index left free in the fit. All point sources with TS>25 

20 in the 6-month source list, lying within the ROI and a surrounding 5°-wide annulus, were modeled 

21 in the fit with single power-law distributions. Different analyses were performed by fitting the 

22 spectra with various models over the whole energy range covered by the LAT above 100 MeV, 

23 or with a PL model over equispaced logarithmic energy bins (where the spectral index was kept 

24 constant and equal to the value fitted over the whole range). In the case of fits with broken 

25 power law (BPL) models, the break energy (EiBreak) bounding the ranges where different photon 
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126 indices (Fi and apply, could not be obtained directly from the fit for most sources because of 

127 convergence problem due the non-smooth character of the BPL function at the break energy. It was 

128 computed from a loglikelihood profile fitting procedure, with statistical uncertainties corresponding 

129 to a dif ference of -2AL= 1 in the loglikelihood (L) with respect to its minimum. We refer the reader 

130 to ref. (|D ' Agostinil |2004| ) regarding limitations with the use of asymmetric uncertainties. 

131 The estimated systematic uncertainty on the flux is 10% at 100 MeV, 5% at 500 MeV and 

132 20% at 10 GeV. The energy resolution is better than 10% over the range of measured E^reak- 



3. Classification 



We employ the co n ventional definition of BL Lac objects outlined in IStocke et al.l (|l99ll ): 
Urry &: Padovani ( 1995 ): Marcha et al. (jl996l ) in which the equivalent width of the strongest opti- 



136 cal emission line is < 5 A and the optical spectrum shows a Ca II H/K break ratio C < 0.4. BLLac 

137 sources were assigned to different subclasses (LSP-BLLacs, ISP-BLLacs and HSP-BLLacs standing 

138 for Low-, Intermediate-, and High-synchrotron peaked BLLacs respectively) according to the posi- 

139 tion of their synchrotron peak, established from radio, optical, UV and X-ray data: fpeafc < 10^^ 

140 Hz for LSP-BLLacs, 10^'^Rz< Up^ak < lO^^Hz for ISP-BLLacs and Up^ak > lO^^Hz for HSP-BLLacs 



141 (jAbdo et al.l l2009a ) . Contemporaneous Swi ft data were used fo r a subset of 46 LB AS sources and 

142 archival ones for the others, as described in lAbdo et al.l ((20091). 



4. Results 
4.1. Photon index distributions 

Although some spectra display significant curvatures, the photon index obtained by fitting 

46 single power-paw models over the whole LAT energy range provides a convenient means to study 

47 the spectral hardness. Fig. [T] displays the distributions of the resulting photon index for the 

48 fo ur different subclass es. The remarkable separation between FRSQs and BLLacs already found 

49 in lAbdo et al.l (j2009d ) is of course still observed for spectra averaged over a 6-month (instead 



50 of 3-month) time span. Likewise, different BLLac subclasses are associated with distinct photon 

51 index distributions in the LAT range. The distributions have (mean, rms)=(2.46, 0.18) for FSRQs, 

52 (2.21, 0.16) for LSP-BLLacs, (2.13, 0. 17) for ISP-BLLacs and (1.86, 0.17) for HSP-BLLacs. For 

53 comparison, the distributions given in lAbdo et al.l (j2009d ) had (mean,rms)= (2.40, 0.17), (1.99, 

54 0.22) for FSRQs and BL Lacs respectively. It must be kept in mind that the 6-month distributions 

55 have been obtained with an improved (more realistic) IRF set leading to a softer measured spectrum 

56 (F higher by ~0.1 unit). Interestingly, the largest difference (increase) in photon index between the 

57 3-month and the 6-month data set is obtained for BL Lacertae (from F=2.24ib0.12 to F=2.54zfc0.07), 



a BLLac intermittently exhibiting broad emission lines characteristic of FSRQs (IVermeulen et al 
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59 ll995l ). It would be interesting to investigate the correlation of this spectral evolution with properties 
60 in other bands. 

The distributions in Fig. [T]are remarkably narrow, with rms =0.16-0.18, i. e. comparable to 

62 the index statistical uncertainties for the faintest LB AS sources. Note that the overlap between 

63 LSP-BLLac and FSRQ photon index distributions is small (LSP-BLLacs being harder), although 

64 both have similar positions of their synchrotron peaks. The gamma-ray photon index is thus a 
.65 distinctive property of a blazar subclass. The low dispersion observed within a subclass strongly 
.66 supports the idea that a very limited number of physical parameters (possibly only one), drive the 
.67 spectrum shape in the GeV energy range. It can also be connected to distinct dominant emission 

68 mechanisms for the different classes. External Compton for the low-ener gy peaked sources and 

69 Synchrotron Self Compton for the high-energy peaked ones, as discussed in lAbdo et al.l (j2009fl ). 



One must keep in mind that the LBAS sample, being significance limited (TS>100 after 3 
71 months of LAT operation) has an intrinsic bias such that faint sources c an more easily be d etected 
if they are hard. This situation is illustrated in Fig. [2] (similar to Fig. 7 in lAbdo et al.ll2009el ) . where 

73 the 6-month average photon index flux is plotted vs the corresponding flux for the four subclasses 

74 along with the approximate LBAS flux limit (for the first 3 months of operation). From this figure, 

75 the source subclass that appears potentially the most affected by this bias is that of H SP-BLLacs, 



76 which are substantially fainter than the other sources. However, as pointed out in lAbdo et al 



77 (|2009a ). more than 60% of known radio- loud HSP-BLLacs are included in the LBAS sample, so the 

78 measured photon index distribution is probably still representative of the full population. 



4.2. Photon index variability 

30 Given the narrowness of the photon index distribution for a given class, the photon index 

31 would not be expected to vary wildly over time for a given source. Two examples of weekly fiux 

32 and photon index light curves are shown in Fig. [3] for the brightest FSRQ (3C 454.3) and LSP- 

83 BLLac (AO 0235+164) in the LBAS. The average photon index is shown as a dashed line in the 

84 corresponding panels. Although flux variations are large (flux variations by a factor >7) for both 

85 sources, the range of photon index values is only about 0.3 wide if one allows for a one-sigma 

86 dispersion. The insets in Fig. [3] display the photon index resulting from an analysis where photons 

87 were sorted in five bins in weekly flux, plotted vs the weekly flux. A weak "harder when brighter" 

88 effect can be seen for both sources. 

To test the constancy of the photon index, the weekly photon indices were fi tted with a constant 



90 model and the corresponding fractional excess variance (IVaughan et al.l l2003l ) was calculated for 

91 the different LBAS sources, keeping only time periods where the sources were detected with a 

92 significance larger than 3 a. Fig. [4] top displays the distributions of normalized values obtained 
.93 for FSRQs (red), BLLacs (blue) and all sources (black). The means of these distributions are close 
94 to 1, as expected for a constant photon index, with no significant differences between FSRQs and 
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195 BLLacs. For illustration, the normalized distribution with 20 degrees of freedom (corresponding 

196 to an average source) expected in the case of a constant photon index is plotted as well. The source 

197 associated with a normalized of 3.5 is PKS 1502+106 for which a clear indic ation of "harder 



198 when brighter" effect has been observed during a bright flare (jAbdo et al.l l2009al ) . Fig. [5] bottom 

199 shows the distributions of fractional excess variance for FSRQs and BLLacs. The means of the 

200 distributions are compatible with (within 1 a); the same applies for all three BLLac subclasses. 

201 The widths of these distri butions are in good agreement with the average statistical uncertainties 

202 (estimated from Eq. 11 in lVaughan et al.ll2003l ). i.e. 8xl0~^ and 1.5x10"^ for FSRQs and BLLacs 

203 respectively. Inspection of distributions of Pearson coefficients of the weekly flux vs photon index 

204 correlation does not reveal any strong trend for any subclass. The average variation of weekly 

205 photon index is plotted vs the relative flux (normalized to the average flux) in Fig. [5] for the 

206 different blazar subclasses. The data from the eight brightest representatives of each subclass have 

207 been considerec|l]. A weak "harder when brighter" effect is apparent in FSRQs, LSP-BLLacs and 

208 ISP-BLLacs, whereas no significant effect is present for HSP-BLLacfl 

209 These observations, aimed at determining the gross spectral features of a source ensemble, do 

210 not exclude fine spectral evolutions over short periods of time, e.g. regarding particular episodes 

211 of fiaring activity, or for particular sources. They do, however, demonstrate that the photon index 

212 in the GeV range changes little with time and within a blazar subclass. 



213 4.3. Spectra of brightest sources 

214 In this section we present the spectra obtained for the eight brightest representatives of the 

215 four subclasses, FSRQs (Fig. E), LSP-BLLacs (Fig. [7]), ISP-BLLacs (Fig. E}, and HSP-BLLacs 

216 (Fig. [9]), ordered according to decreasing average flux. Upper limits are shown for bins associated 

217 with a TS lower than 9 (significance lower than 3 a) or with a number of source photons (predicted 

218 by the model) lower than 3. The brightest FSRQ, 3C 454.3 with an average flux (E>100 MeV) o f 



219 2x 10~^ ph cm~^s~^, exhibits a pronounced break around 2 GeV, as reported in lAbdo et al.l (|2009dl ). 

220 Indications for breaks between 1 and 10 GeV are observed for essentially all of these FSRQ sources. 

221 This behavior is conflrmed by comparing (Fig. [TT]) the flux (E>2 GeV) extrapolated from the 

222 spectral distribution in the range 100 MeV <E< 2 GeV, fltted with a power-law function, with the 

223 actual measured flux (obtained via a power-law flt in the E>2 GeV energy range). For all sources 

224 considered, the measured flux is lower than the extrapolated flux by more than 30%. The spectral 

225 properties of these sources are summarized in Table 1, which gives also the difference in loglikelihood 

226 between the BPL and PL flts, the break energy in the source rest frame, Brpn.k=^Rrpn.k>i(^+z), 

227 and the source gamma-ray luminosity (computed as in lGhisellini et allbood ). E' Break lies between 



condition on the significance > la, i.e. less restrictive than above, has been imposed in this analysis to minimize 
the effect of the instrumental bias against soft and faint states illustrated in Fig. (2] 

^The weak tendency toward softer spectra at high flux in HSP-BLLacs does not appear to be significant. 
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228 1.6 and 10 GeV for all sources listed in Table 1. No significant correlation is found between Break 

229 and the gamma-ray luminosity (Fig. [TO]) . 



230 The presence of a break is also clear for two of the brightest LSP-BLLacs (associated with 

231 the highest luminosities, Table 2) in the LBAS sample, AO 0235+164 and PKS 0537-441 (Fig. W), 

232 while it is less apparent for the fainter ones. From the inspection of these spectra, it can already 

233 be noted that the onset of the break, when present, seems to be located at higher energy than in 

234 the case of FSRQs while E'^^eafc lies in about the same range. 

235 Some ISP-BLLacs (Fig. Ej) present clear signs of breaks (e.g. S5 0716+71, S2 0109+22), while 

236 the rest, having in some cases very hard spectra, are compatible with power-law distributions up to 

237 several tens of GeV. Finally, no bright HSP-BLLac (Fig. [9]) shows any evidence for a break in the 

238 LAT energy range. This simple observation, which could be naively expected for sources many of 

239 which are detected at TeV energies, definitely rules out an instrumental effect as the origin of the 

240 break found for lower-e nergy peaked sources. The indication for a moderate break in PKS 2155-304 



241 (jAharonian et al.ll2009l ) observed over a short time period (11 days) does not persist over a 6- month 



242 integration time. 

243 For a small number of sources exhibiting the break, a few photons compatible with the source 

244 location are detected at high energy, at variance with the decreasing trend. Most of these have 

245 been cut off by the condition on the minimum number of photons per bin (3). More statistics 

246 will be required to determine whether these photons do arise from these sources or just represent 

247 background fluctuations. 

248 Fig. [TT] illustrates the general trend for the four subclasses. The trend observed in the presence 

249 of a spectral break (or curvature) parallels that observed in the photon index for the four different 

250 classes. 



4.4. Spectra of special sources 



252 For 19 of the 22 brightest LBAS FSRQs, a likelihood ratio test (LRT, iMattox et al.l \19M) 

253 rejects the hypothesis that the spectrum is a PL (null hypothesis) against the one that the spectrum 

254 is a BPL, at a confidence level greater than 97%. The four top panels of Fig. [12] show spectra of 

255 representative sources where the break is clear. The spectral properties of these sources (PL and 

256 BPL fit results, difference in logLikelihood between the two fits) are reported in Table 1 as well. 

257 The two bottom left panels correspond to two of the three sources having the confidence level less 

258 than 97% (the other being PKS 2022-07, whose spectrum is given in Fig. [6|). For these two sources, 

259 a break located around 10 GeV cannot be excluded. 

260 The last two panels of Fig. [12] present the energy spectra of sources of particular interest, 

261 Mrk501 and lES 0502+675. Mrk501, the archetypal example of an extreme HSP-BLLac, is hard 

262 (r=1.75±0.06) and does not exhibit any sign of curvature, in keeping with the behavior of the 
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263 other HSP-BLLacs. The other one corresponds to lES 0502+675, an HSP-BLLac which exhibits 

264 an unusual concave energy spectrum (ri=2.68±0.18,r2 =1.47+0.10, Eiyr-eak= 1.4+0.6 GeV). The 

265 LRT indicates that a BPL model is favored against a PL with a significance level of 2 x 10~^. The 

266 very hard spectrum above 1.4 GeV would make this source a prime target for TeV observations 

267 although the redshift is fairly large (0.416). 

268 4.5. Detailed analysis of 3C 454.3 and AO 0235+164 energy spectra 

269 The data accumulated over 6 months enable us to discriminate between different spectral 

270 models for the two brightest sources with spectra exhibiting strong departure from a pure PL, 

271 namely 3C 454.3 and AO 0235+164. Fig. [T3l shows the results of fits with different models: PL 

272 (thin lines), BPL(thick solid) and log-parabola (dashed) are compared with the data. Surpris- 

273 ingly, a broken power-law model is favored as the best fit for both sources. Despite the good 

274 statistics, no curvature is apparent in the energy range below the break. The fitted parameters 

275 are: F[E>100 MeV]=1.97+0.03x 10"^ ph cm-^s"!, ri=2. 39+0.02, r2=3. 42+0.11, Ebreak= 2.5+0.3 

276 GeV for 3C 454.3 andF[E>100 MeV]=0.60+0.02x 10"^ ph cm^^g-i, ri=2. 05+0.02, r2=2. 95+0.16, 

277 Efereafc=4.5l};oGeV for AO 0235+164. While the break energy is somewhat larger for AO 0235+164 

278 than for 3C 454.3, the photon index change (r2 - Fi) is about the same (0.90+0.16 vs 1.03+0.11). 

279 The similarity of the break feature for two sources belonging to different subclasses, FSRQ and LSP- 

280 BLLac, with different line strengths, seems to rule out any absorption effect in the LSP-BLLacs. 

281 4.6. Apparent curvature due to varying spectral hardness 

282 An energy spectrum with time varying hardness may exhibit an apparent curvature when in- 

283 tegrated over an extended period of time. To assess the magnitude of this effect, both analytical 

284 estimates and simulations assuming pure power-law distributions with flux and photon index cor- 

285 responding to those actually measured over weekly time bins have been performed. Fig. [3] shows 

286 the 3C 454.3 and AO 0235+164 weekly light curves and corresponding photon index. With these 

287 input data, the calculated spectra (assuming constant exposure) are not found to exhibit significant 

288 curvatures (Fig. [HI) . Fig. [15] compares the simulated photon count distribution obtained within 

289 the 90% containment radius around the source (blue) with the data (black). The effect of spectral 

290 hardness varying with time clearly cannot alone account for the observed features. 

291 5. Discussion 

292 The trend in the observed gamma-ray photon index reported here confirms that reported earlier 

293 using three months of data: FSRQs, with a gamma-ray photon index greater than 2, are softer than 

294 BLLacs, indicating that the peak of the high-energy component in FSRQS is always lower than 
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295 100 MeV. For BLLacs, the gamma-ray photon index correlates with the different BLLac subclasses, 

296 which themselves are defined by the position of the synchrotron peak. The measured photon index 

297 shifts from T >2 to F <2, indicating that the peak energy of the high-energy component of the 

298 spectral energy distribution (SED) sweeps across the Fermi energy range from LSP-BLLacs to 

299 HSP-BLLacs. 

300 The photon index being related to the shape of the emitting electron energy distributions, 

301 different regions of the electron distributions (from the high- to the low-energy ends for FSRQs/LSP- 

302 BLLacs and HSP-BLLacs respectively) are probed in the LAT range. One would thus expect 

303 different spectral variability patterns in the LAT energy range for different blazar classes. 

304 Using weekly light curves obtained over the first 6 months of LAT operation, this expectation 

305 does not seem to be corroborated by the data. The gamma-ray photon index appears remarkably 

306 stable with time, irrespective of the blazar class. This feature is in line with the observed narrowness 

307 of the index distribution for a class: a strongly varying photon index for a source would inevitably 

308 lead to a broad class distribution as different sources would be "caught" in different states. This 

309 apparent constancy of the photon index for all blazar classes may appear surprising, since for 

310 different classes, different electron energies (potentially associated with different cooling timescales, 

311 tcooi) emit gamma-rays in the LAT energy range. However, estimates of the cooling time show that 

312 within one week, all pairs with Lorentz factors 7 > 5/6iou' have cooled down (where 610 is the 

313 Doppler factor 6 /{I + z) divided by 10, and u' is the jet frame energy density in erg cm~^ of the 

314 ambient magnetic field, or photon field provided scattering occurs in the Thomson regime). For 

315 sufficiently large u' the LAT energy range falls in the complete cooling regime. Spectral hysteresis 

316 may then be expected on time scales shorter than weekly if the duration of electron injection is 

317 sufficiently limited <^ tcooi- Flux variations on weekly time scales may then reflect a varying injected 

318 energy content into the emission region from one week to another. Alternatively, continuous particle 

319 injection on at least weeks time scale could stabilize the spectral index. In the case of lower u', 

320 constant spectral indices can still be expected during fiux variations in the LAT energy range if 

321 continuous particle injection is maintained for a time range significantly lo nger than the cooling 



322 time of the gamma-ray emitting particles (e.g. iKirk &: Mastichiadid Il999l ). Flux changes then 

323 constrain the duration of particle injection. As more data are accumulated, possible exceptions to 

324 the observed index stability may appear. However, we can safely claim that a ver y soft spectrum 



325 for a n HSP-BLLac, hke that reported from EGRET for PKS 2155-304, F=2.35 (jHartman et al 



19991). or a very hard one for a FSRQ, like the one reported from AGILE for 30 454.3, r=1.7 



327 (jVercellone et al.ll2009l ). represent very rare occurrences. 



328 All blazar spectra measured by EGRET were represented with pure PL. Thanks to its improved 

329 sensitivity, Fermi has revealed that the spectra of some low-energy peaked blazars display strong 

330 departure from a pure-PL behavior, with a BPL function as the best model. This feature being 

331 present for essentially all FSRQs where it can be detected with sufficient significance, it is likely to be 

332 a general character of this class. The three brightest LSP-BLLacs, in which this effect is also clearly 

333 seen, have higher luminosity (L-y >10'^^ erg s~^) than the rest of the LBAS LSP-BLLacs. Two of 
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334 them are known to exhibit broad emission hnes in low emission states, s o could be FSRQs w hose 

335 lines are hidden by non-thermal emission in active states. As discussed in lAbdo et al.l (j2009dl ). the 

336 difference in photon index for most sources is significantly larger than 0.5 expected for an incomplete 

337 cooling effect. An absorption effect seems to be ruled out as well, since to produce a break in the 1- 

338 1 GeV, the photon field should have an energy peaking in the 0.05-0.5 KeV range, which excludes 

339 the broad-line region peaking in the UV. This feature most likely reflects the energy distribution 

340 of the emitting electrons. For the low-energy peaked sources where it is seen, T is greater than 

341 2, i.e. the LAT range corresponds to "the falling edge" of the IC hump, where the highest energy 

342 electrons contribute. This feature could indi cate a cutoff in that distributio n, possibly related to 



343 limitations in the acceleration process (e.g. iDrurv 



1991 



Webb et al.l 1 19841 ) . Comparison of the 



344 maximum electron energies determined by the SED for PKS 2155-304, an HSP-BLLac where a 

345 one-zone synchrotron/SSC model was applied to derive the magnetic field and Doppler factor, 

346 showed that such an interpretation requires that the acceleration rate be approximately 3 orders 

347 of magnitude sma ller than the maximum acceleration rate determined by the Larmor timescale 

348 (jFinke et al.l l2008l ) . A Klein-Nishina break, expected in the context of a dominating External 

349 Compton emission process where the electrons upsc atter BLR photons, is predict ed to set in around 

350 I'^i^z) GeV, where F is the blob bulk Lorentz factor ( Ghisellini &: Tavecchidl2009l ) . i. e. significantly 

351 higher than the energies found here. 

352 Irrespective of its origin, this feature, common for FSRQs and some LSP-BLLacs, has impor- 

353 tant practical consequences. First, it surely complicates the assessment of EBL attenuation effects 

354 using FSRQs and LSP-BLLacs, as fewer photons are detected in the >10 GeV energy range. Sec- 

355 ond, as the low-energy peaked blazars are likely to represent the bulk of the blazar population, this 

356 break must manifest itself in (and be considered when evaluating) the contribution of unresolved 

357 blazars to the extragalactic diffuse 7-ray background. Finally, this effect must be considered when 

358 estimating the detectability of a source in the TeV range. 

359 The concave shaped SED measured for the X-ray selected object lES 0502+675 (z=0.416) 

360 opens interesting questions. It could potentially be a spurious feature resulting from the spatial 

361 confusion of a hard source with a soft one. No evidence for a second source has been found using 

362 the gttsmap tool of the standard ScienceTools package. The closest CRATES source is 1.5° away. 

363 If confirmed, this peculiar spectral shape indicates either two components (e.g. Synchrotron-Self 

364 Compton and External Compton in the context of leptonic models) in the high-energy SED (how- 

365 ever not expected in HSP-BLLacs in the framework of one- zone leptonic models), or the turnover 

366 from the synchrotr on to the high-ene rgy component, although the contemporaneous spectral en- 

367 ergy distributions (jAbdo et al.ll2009fl ) obtained during that period does not seem to support that 

368 interpretation. This object should be a prime target for TeV instruments, opening interesting 

369 perspectives for studies on EBL absorption at TeV energies given its redshift. 

370 A similar behavior, although less significant, has also been found for lES 1959+650 (Fig. [8]) 

371 and PC 1246+586, which is not included in the LBAS sample but is present in the LAT Bright 
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372 Source (jAbdo et al.ll2009d ). The spectra of other sources (e.g. W Comae) exhibit a wavy 



373 shape possibly indicative of multiple components as mentioned above. More detailed analysis will 

374 be necessary to resolve this issue. 



375 6. Conclusion 

376 The spectral properties of the LBAS blazars in the gamma-ray band, as determined over 

377 the first 6 months of LAT operation, have been presented. The average photon index of LBAS 

378 blazars are found to be r=2.46 for FSRQs, r= 2.21 for LSP-BLLacs, r=2.13 for ISP-BLLacs 

379 and r=1.86 for HSP-BLLacs, with an rms of 0.16-0.18. Spectral breaks have been observed to 
3B0 be common features in FSRQs (the break energy ranging from 1 to 10 GeV in the source frame 
3B1 for the brightest sources), and present also in some bright LSP-BLLacs. The different spectral 

382 features reported here represent challenges for theoretical models aiming at describing the blazar 

383 phenomenon. Althou gh the fairly strong correlation between photon index and blazar class fits 



384 well within pictures (iGhisellini et al.l Il998l : iBottcher &: Dermeij |2002| ) where the cooling due to 



385 strong ambient radiation fields (manifesting themselves via the presence of emission lines) limits 

386 the acceleration of particles at high energy, the near constancy of this photon index with time and 

387 flux variation provides new constraints on the emitting particle dynamics. Moreover, the fact that 

388 spectra for most FSRQs and some LSP-BLLacs are best modeled by a BPL with a break in the 

389 1-10 GeV range is quite unexpected, the break representing a distinctive feature of these sources. 

390 The LAT has already revealed novel aspects of gamma-ray blazars and will help refine the new 

391 picture that progressively emerges as more data (both in the gamma-ray and other bands) are 

392 accumulated. 
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Fig. 1. — Gamma-ray photon index distributions for the four blazar subclasses. 
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Fig. 2. — Photon index T vs Flux (E > 100 MeV) for the LB AS sources considered here. Open 
circles, red: FSRQs; solid symbols: BL Lacs (cyan circles: LSP-BLLacs, green triangles: ISP- 
BLLacs, blue squares: HSP-BLLacs). The solid curve represents the TS = 100 limit for a 3 month 
period (i.e. the defining condition of the LBAS sample) estimated for (l,b)=(80°, 40°). 
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Fig. 3. — Measured weekly fluxes and photon index for 3C 454.3 and AO 0235+164 used as input of 
the simulations described in the text. The insets show the photon index resulting from an analysis 
where photons were sorted in five weekly-flux bins plotted vs the weekly flux. 
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Normalized excess variance 

Fig. 4. — Top: Normalized distributions of the weekly photon indices for FSRQs (red), BLLacs 
(blue) and all sources (black). Bottom: Normalized excess variance distributions of the weekly 
photon indices for FSRQs (red) and BLLacs (blue). 
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Fig. 5. — Average variation of photon index, F, vs relative flux (normalized to the average flux), 
for the different blazar subclasses. The data are integrated over a weekly timescale, and only the 
eight brightest representatives of each subclass have been considered. 
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Fig. 6. — Gamma-ray spectra of the eight brightest FSRQs in the LBAS sample obtained for 
equispaced logarithmic bins (dots), together with residues with respect to the fitted power-law 
model (solid lines in upper panels). 
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Fig. 9. — Same as Fig. [H for the eight brightest HSP-BLLacs in the LBAS sample. 
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Fig. 10. — Break energy vs gamma-ray luminosity for the FSRQs (closed red symbols) and LSP- 
BLLacs (open blue symbols) listed in Table 1 and 2 respectively. 
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Fig. 11. — Measured vs extrapolated flux above 2 GeV for the eight brightest sources of each blazar 
subclass. The distance of the points from the diagonal is indicative of the presence of a spectral 
break. 
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Fig. 12. — Same as Fig. [U for eight particular sources in the LBAS sample. 
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Fig. 13. — Upper panels: Energy spectra of 3C 454.3 (left) and AO 0235+164 (right) compared 
with fit results obtained with different models: PL (thin solid), BPL (thick solid) and logparabola 
(dashed). 
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Fig. 14. — Analytical energy spectra (solid curves) of 3C 454.3 and AO 0235+164 resulting from 
summing power-law distributions with parameters (flux, photon index) as measured in weekly bins 
(Fig. 2]). The dashed lines represent PL fits. 




Fig. 15. — Count distributions within the 90% containment radius of simulated (blue, thin) and 
real (black, thick) data for 3C 454.3 and AO 0235+164. The solid curve correspond to the total 
contribution of galactic and isotropic diffuse backgrounds. 



Table 1. Spectral properties of selected FSRQs 



Name 


1 


b 








F 




AL 








Fi 






F2 






AF 






z 


E'B^ofc (GeV) 


Luminosity^ 


3C 464.3 


86.12 


-38.1 


2, 


.a53±a. 


.02 


2, 


.47±0, 


.01 


-64.7 


1 


.994±0, 


.029 
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.40±0, 


.01 


3 


.61±0. 


12 


1. 
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,12 


2 5+" =* 
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4g+0.55 
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44.1 


PKS 1502+106 
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1, 
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1 
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3, 


.06±0. 


12 


0. 


,89±0. 


,12 
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1.839 


"■'-0.4 
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57.06 


0, 
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.01 


2, 
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0, 
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2, 
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.03 


2, 


,60±0, 


,06 


0. 


.25±0. 


.06 




0.536 




5.0 


PKS 1610-08 


361.2 


40.13 


0, 


.739±0, 


.02 


2. 


.47±0. 


.02 


-7.13 


0, 


.717±0, 


.042 


2, 


.42±0, 


.06 


3 


.08±0. 


26 


0. 


.66±0. 


,26 




0.36 


„ n+0.9 


1.60 


3C 273 


289.9 


64.36 


0, 


.682±0, 


.02 


2. 


.73±0, 


.03 


-9.12 


0, 


.669±0, 


.023 


2, 


.68±0, 


.03 


3 


.66±0. 


,28 


0. 


,97±0. 


,28 




0.158 


1 n+0.55 
^■"-0.55 


0.19 


PKS 0454-234 


223.7 


-34.9 


0, 


,632±0, 


,01 


2, 


,19±0, 


.01 


-23.6 


0, 


,604±0, 


.016 


2, 


.11±0, 


.02 


3, 


.28±0. 


21 


1, 


,16±0, 


,21 




1.003 


9 95+2'l 
».»o_Q 25 


22.1 


PKS 2022-07 


36.89 


-24.3 


0, 


,439±0, 


,01 


2 


,38+0, 


,03 


-4.54 


0, 


,420+0, 


,018 


2 


,32+0, 


,03 


2 


,84+0. 


17 


0, 


,52+0, 


18 


^' '"'-0-18 


1.388 


., + 1-8 
''•*'-0-45 


31.1 


TXS 1520+319 


50.14 


57.04 





,381±0 


,01 


2 


,52+0, 


,03 


-3.95 


0, 


,364+0, 


,016 


2 


,45+0, 


,04 


2 


,90+0. 


15 


0, 


,45+0, 


16 


-'■°-0.4 


1.487 


4 0+^-" 
*'"-0-l 


35.4 


4C +38.41 


61.08 


42.40 





.229+0 


,01 


2 


,50+0, 


,04 


-10.1 


0, 


,211+0, 


,014 


2 


,36+0, 


,05 


4, 


,()6+(). 


50 


1 


70+0, 


,50 


, ., + 0-55 
^■•^-0-5 


1.814 


6 4+1-'' 
"■*-0-15 


39.3 


PKS 1908-201 


16.89 


-13.1 





,160+0 


,01 


2 


,42+0 


,06 


-9.31 


0, 


,143+0 


,018 


2 


,09+0 


,13 


3, 


,18+0. 


22 


1 


,09+0, 


,25 


'^■^-0.25 


1.119 


^■"'^-0.55 


7.75 


PKS 0528+134 


191.3 


-11.0 


0, 


.309+0 


,02 


2 


,72+0, 


,05 


-7.42 


0, 


,268+0 


,020 


2 


,52 + 


,07 


3, 


,87+0. 


47 


1, 


,34+0, 


,48 


J g+0.45 


2.07 




71.4 


RGB J0920+446 


176.7 


44.82 


0, 


.247±0, 


,01 


2, 


,22±0, 


,03 


-7.61 


0, 


.235+0, 


.012 


2, 


.16+0, 


.03 


4, 


.86+0. 


,84 


2, 


,70+0, 


,84 




2.19 


„ „ + 11.0 
■^■■^-1.5 


67.2 



»10-'' ph[B>100 MeVJcm-^s-l 
'■lO^'' erg s-1 



Table 2. Spectral properties of selected BLLacs 



Name 


1 


b 




r 


AI, 


FluxBPi" 


Pi 


r2 


AP 


BBreofe (GeV) 


z 


E/Br-eak (GeV) 


Luminosity^ 


LSP-BLLacs 
AO 0235+164 
PKS 0537-441 
PKS 0426-380 


156.7 
250.0 
240.7 


-39.0 
-31.0 
-43.6 


0.630±0.01 
0.400±0.01 
0.274±0.01 


2.12±0.01 
2.28±0.02 
2.18±0.03 


-20.8 
-8.51 
-7.17 


().599±().()16 
().38()±().()15 
0.255±0.012 


2.()4±().()2 
2.20±().()3 
2.10±0.03 


2.80±0.12 
3.08±0.23 
3.36±0.45 


0.75±0.13 
0.87±0.24 
1.26±0.45 


4 5+1-^ 
^■''-o.s 

o , + 6.0 
°-^-0.4 


0.94 
0.892 
1.112 


K S+3-0 

'■^-1.5 
^+13.0 


19.6 
10.1 
12.3 


HSP-BLLac 
lES 0B02-H675 


143.7 


15.89 


0.019±0.00 


1.70±0.14 


-8.40 


0.064±0.015 


2.68±a.l8 


1.47±0.10 


-i.2±a.2i 


1 4+°-'' 


0.416 


2 0+^° 
"'■"-0.7 


0.12 



^lO-'' ph[E>100 MeV]cm-^s-l 
'•10*'' erg 



